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A solvable vertex model in ferromagnetic regime gives rise to a soliton cellular 
automaton which is a discrete dynamical system in which site variables take 
on values in a finite set. We study the scattering of a class of soliton cellular 
automata associated with the [/ g (.Dn )-perfect crystal B 2,s . We calculate the 

2,s o d2,1 



combinatorial R matrix for all elements of B ' s <g) B > . In particular, we show 



that the scattering rule for our soliton cellular automaton can be identified 
with the combinatorial R matrix for U q (A^) © ^(-D^j^-crystals. 



I. INTRODUCTION 

A cellular automaton is a dynamical system in which points in the one- 
dimensional space lattice are assigned discrete values in a finite set which evolve 
according to a deterministic rule. Soliton cellular automata (SCA) are a kind 
of cellular automata which possess stable configurations analogous to solitons 
in integrable partial differential equations. Solitons move with constant velocity 
proportional to the length if there is no collision. After collision their lengths 
are preserved but the phases are shifted. Moreover, SCAs have many conserved 
quantities. The simplest example of a SCA is the Takahashi-Satsuma's automaton 



3 If ] where the site variables take on two values {0, 1}. This SCA gives rise to a 
nonlinear dynamical system. Subsequently, this was generalized to other SCAs 
where the site variables take on more than two values ([30|. [32|. [33fJ ) . It was then 
found that these systems can be described by perfect crystals [15| for the quantum 
affine algebra U' q {A^) (8, @). In fact, it was shown in |4] that the phase shift of 

the U' q {A^) SCA was related to the energy function of U' q (An^) perfect crystals. 
In jof, a class of SCAs associated with the perfect crystals for the non-exceptional 
quantum affine Lie algebras given in [3| was constructed. In (35[ (resp. [26]) the 
SCA associated with the U' q (D^) (resp. U' q (G 2 )) perfect crystals given in [l7[ 
(resp. j25|) have been constructed. In this formulation the time evolution operator 
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is given by the row-to- row transfer matrix of integrable vertex models at q — 0, 
represented as a product of the combinatorial R matrix of perfect crystals. Using 
this approach, the scattering rules for the SCA associated with the perfect crystals 
given in [14] were determined in [Io|]. Among these SCAs, the SCA associated 
with the quantum affine algebra U'(Dn ) is of fundamental importance, because 
the other ones can be imbedded into this one [20| . 

It is known that the perfect crystals for quantum affine Lie algebras are crystals 
for certain Kirillov-Reshetikhin (KR) modules [l9[. Fourier, Okado and Schilling 
(0] ,0) have proved the existence of perfect crystals for quantum affine algebras 
of classical types conjectured in (jsj, j9[). In particular, in j2[ explicit descriptions 
of KR-crystals B r,s for any Dynkin node r ^ and any positive integer s are given 
for each quantum affine algebras of classical types. In [3j it has been shown that 
if s is a multiple of t r = max{l, 2/ (a r , a r )}, then the KR-crystal B r,s is perfect of 
level s/t r . So far all SCAs constructed are associated with the perfect crystal B 1,a 
with the exception of the one in j24j where the SCA associated with the L^(-Di^) 
perfect crystal B n ' s has been constructed. It is conjectured that the scattering rule 
for the SCA associated with a perfect crystal B r ' s will be given by the R matrix for 
the crystal of the quantum affine algebra U'Jg'^ 1 ') © U' q (g"^), where g' and g" are 
the simple Lie algebras associated with the two connected components obtained 
after removing the Dynkin node r from the finite Dynkin diagram. 

In this paper, using the realizations of B 2,s given in [2} (see also [28]) we compute 
the combinatorial R matrix B 2 ' s © B 2 ' 1 ^B 2 ' 1 © B 2yS for U' q (D^') : then use it to find 
the scattering rule of the SCA associated with B 2,s . We see that solitons of length 
s are parameterized by the U q (A^) © U q (D^}_ 2 )-cryst8l B x,s x B 1,s . Furthermore, 
when two solitons collide, the scattering rule is described by the combinatorial R 
matrix B 1,S1 © B 1,S2S ^B 1 ' S2 © B 1,S1 , si > S2- Thus we have shown that the above 
conjecture holds in this case. 



II. PRELIMINARIES 



The infinite dimensional Kac-Moody algebra g = Dn^ is determined by the 
following Dynkin diagram: 

1 n-1 



2 3 



n — 2 n 



We let P = span z {A , A 1; . . . , A n } be the weight lattice, {a , cei, ■ ■ ■ , a n } be the 
set of simple roots, and {h , hi, . . . , h n } be the set of simple coroots. 
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A. Crystals 

In this section we give the basic definitions regarding crystals. Our exposition 
follows [Hi. Let / = {0, 1,2, ...n}. 

Definition 1. A crystal associated with U q (o) is a set B together with maps wt : 
B — > P,e~i, fi : B — > B U {0}, and Ei,(fi : B — > Z U {— oo} ; for i G / satisfying the 
following properties: 

1. (fi(b) = Ei(b) + (hi, wt(h)) for all i E I, 

2. wt{eib) = wt{b) + c*i if eft G B, 

3. wt(fib) = wt(b) - ^ iffibeB, 

4. Ei(eib) = Ei(b) - l,<fi(eib) = (fi(b) + 1 if e { b G B, 

5. EiCfib) = Ei{b) + l,ipi{fib) = <pi(b) - 1 if fib G B, 

6. fib = b' if and only if b = e~ib' for b,b' G B,i G /, 

7. if (fi(b) = —00 for b G B, then e~ih = fib = 0. 

A crystal B can be regarded as a colored, oriented graph by defining 

^ fib = b'. 

Definition 2. The tensor product Bi®B 2 of crystals B\ and B 2 is the set B\ x B 2 

together with the following maps: 

1. wt(b! <g> b 2 ) = wtib^ + wt{b 2 ), 

2. Eiipi ® b 2 ) = max(£j(6i), £j(& 2 ) - (h h wtipx))), 

3. <fi(bi <g> b 2 ) = max(^(6 2 ), fiih) + {hi, wt(b 2 ))), 



4. Ci(bi 



5. fi{h 



Cib x ® b 2 , if<fi(bi) > £i(b 2 ), 

h <S> eib 2 , if<fi(h) < £i{b 2 ), 

fih ®b 2 , ififiih) > Ei(b 2 ), 

h <g> fih, if<fi(h) < Ei(b 2 ), 



where we write bi <S> b 2 for {pi, b 2 ) G B\ x B 2 , and understand b\ ® = ® b 2 = 0. 

B\ £g> B 2 is a crystal, as can easily be shown. 

Definition 3. Let B\ and B 2 be U q {o)- crystals. A crystal isomorphism is a bijec- 
tive map \P : B v U {0} -)■ B 2 U {0} such that 
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1. (0) = 0, 

2. if be B 1 and *(6) G B 2 , then wt(V(b)) = wt{b),ei{^{b)) = e^b), ipi(^(b)) = 
(fi(b) for all i £ I, 

3. ifb,b' G S 1 ,*(6),*(6') G 5 2 and fjb = b', then f^{b) = *(6 / ) and #(&) = 
e^(o') /or a// i G J. 



B. -crystal 5 M 

Let B 1 ' 8 := {6 = (x 1 , x 2 , . . . , x n , x n , x n _i, . . . , x x ) G Z| n |s(6) := XT=i x * + 



SiLi x i = s ? x n = or x n = 0} and define 
e b 



e n b 



lib 



fob 

fnb 

fib 



(xi, x 2 - 1, • • • , x 2 , xi + 1) if x 2 > x 2 , 
(x 1 — l,x 2 ,...,x 2 + l, xi) if x 2 < x 2 , 

(xi, . . . , x n —\ -\- 1, x n , x n 1, . . . , X\) if x n 0, x n > 0, 

(xi, . . . , Xi -\- 1, Xj_|_i 1, . . . , xi) if Xj-fi i> Xj_|_x, _^ ^ ^ ^ _^ 
(xi, . . . , x i+ i + 1, Xi - 1, . . . , xi) if 

(xi, x 2 + 1, . . . , x 2 , xi - 1) if x 2 > x 2 , 
(xi + 1, x 2 , . . . , x 2 - 1, xi) if x 2 < x 2 , 

(x^, . . . , x n 1, x n; x n —\ ~\~ 1, . . . , x^) if x n ^> 0, x n 0, 
(x^, • • • •} x n —\ 1, x n; x n ~)~ 1, . . . , X\) if x n 0, x n ^ ; 

(xi, . . . , Xi — 1, Xj + i + 1, . . . , Xi) if Xj +1 > x i+1 , j ^ ^ _ ^ 
(xi, . . . ,x m - l,Xi + 1, . . . ,xi) if x i+1 < x i+1 . ~ ~ 



If Xi < or Xi < in b' = ej(6) or /j(6) then b' is understood to be 0. 



n-2 



wt(6) = (xi - xi + x 2 - x 2 )A + y^(xj - Xi + Xj+i - x i+ i)Ai 

i=i 

"I - (-^n— 1 -^n— 1 "I - %n •En)A-n—l 
H~ (^n— 1 *^n— 1 "I - 2'ra)A. n , 

V? (o) = xi + (x 2 - x 2 ) + , e (b) = x 1 + (x 2 - x 2 ) + , 
V?j(6) = Xi + (xj+i - x i+ i)+ for 1 < i < n - 2, 
£j(6) = Xi + (x i+ i - x i+ i) + for 1 < i < n - 2, 
y?„_i(o) = x n _i + x n , £ n -i(&) = x n -i + x n , 

where (n) + := max(n, 0). Then we have the following: 
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Theorem 1 ([l3],[lH). The maps e~i, fi,Ei,f>i,wt define a U'(DfP) -crystal struc- 
ture on B 1,s . 

We associate the element (xi,x 2 , ■ ■ ■ ,x~i) G -B 1 '* with the tableau: 



1 


1 




1 


2 


2 




2 




1 


1 




1 














XI 








X2 













C. Perfect crystal B~> s for Dn 



In this section we review the perfect crystals B 2,s , s > 1 corresponding to the 
2- node of the Dynkin diagram of D„ ■ The existence of crystal bases B 2 ' s ,s > 1 
was proven in [23] and the combinatorial realization was given in [2[ (see also 28 

0)- 

For D { n\ define the alphabet B = {1,2, 
the following partial ordering on B: 



n — 1, n, n, n — 1, 



2,1}. Define 



l<2<---<n — 1< _ <n 

n 



1 < •■■ < 2 < 1. 



For D n , define the set: 



B(kA 2 



T 



Ti,i 


Tl,2 




Tx, k 


T2.1 


T^2,2 







rp ^ rp 



l,2,j = l,2,...fc-l 



T\,j < T 2>j or Tij 
3 = 1,2,. ..k, 

no column of the form 



occurs, 



no configuration of the form 
occurs, and 



a 


a 


* 


a 



or 



a 


* 


a 


a 



no configuration of the form 
or 



?2—l 


n 


n 


n — 1 



n — 1 


n 


n 


n — 1 



occurs. 



The set B(kA 2 ) becomes a C/'(D n )-crystal by considering it as a tensor product 
(see Definition [2]) under the following reading: 



Ti,i 


Tl,2 






T2,l 


T 2i2 




T 2) k 
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corresponds to 

T ltk <g> T 2jfc ® ri )fc _i ® T 2ife _! ® • ■ • (8) Ti,i (8) T 2jl . 
The maps ej, fi,Ei,tpi,i G / and wt are given for an individual letter b E B by: 





' i, if b = i + 1, 








ei{b) — < 


z + 1, it b — %, i = 1, 
k 0, otherwise. 

n, if o = n — 1, 


2, ... ,71 


- -L> 




e n (b) = < 


n — 1, if 6 = n, 
, 0, otherwise. 

f i + 1, if 6 = z, 






(2) 


h(b) = < 


i, if 6 = i + 1, i — 1, 
k 0, otherwise. 


2, .... n 


- 1, 


(3) 




ri — 1, if 6 = n, 








/»(&) = < 


n, if 6 = n — 1, 
k 0, otherwise. 






(4) 


1 

£i(6) = < 


1 . if b — i 4- 1. i. 
v 0, otherwise. 


...,7l- 


1, 


(5) 


e n {b) = < 


1, if 6 = n — 1, n, 
v 0, otherwise. 






(6) 


^(6) = < 


' 1, if 6 = + 1, f _ 1 2 
v 0, otherwise. 


...,n - 


1, 


(7) 


<Pn(b) = < 


1, if b = n — 1, n, 
^ 0, otherwise. 

'A!, if 6 = 


1, 




(8) 


wt(6) = < 


Ai-Aj-i, if 6 = 
, A n _i + A n - A n _ 2 , if & = 


2,3,... 
n — 1, 


71 — 2, 71, 


(9) 


wt(6) = - 


-wt(6), in any other case. 









We then have the following: 

Theorem 2 ((27j],|2]). B 2 ' s ,s > is a D$P -crystal, and, forgetting the 0-arrows, 
is isomorphic to the D n -crystal Q s k=Q B(kA2). 

The 0-action in B 2,s was described in [28[ and later in [2[ for arbitrary B r,s . It 
is given in terms of the .D^-crystal automorphism a, induced by the symmetry 
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between the and 1 nodes in the Dynkin diagram. We give a method for computing 
er, and use it to give expressions for e , f ,e ,(p . 

Consider the restriction of B(kA 2 ),0 < k < s from D n to -D„_i by delet- 
ing the 1 arrows. The D n _i branching components are the irreducible com- 
ponents of the restricted crystal. These may be partially ordered by setting 
B < C if there exists a 1-arrow connecting some element of B to some ele- 
ment of C in B(kA 2 ). There is also a rank function relative to s and k asso- 
ciated with each branching component B given by setting ik(B) = #{C\B < 
C in some path from B to the highest branching component in B(kA 2 )} + s — k + 
1. We define the function i) : B(jA 2 ) U {0} ->■ B(kA 2 ) U {0}, < j, k < s by defin- 
ing Lj(T) to be the element corresponding to T in the corresponding branching 
component of B(kA 2 ) whose rank is equal to that of B, if it exists, and if no corre- 
sponding branching component exists, and t J k (0) = 0. Then, ^ +1 (T), < k < s — 1 
may be given explicitly. Let T' be the 2 x k + 1 tableau defined by the following: 
if k ^ 0, then 




and, for 2 < i < k 




2,i - 




T 2> i, if T 2)i 7^ Ti^-i 
n, if T 2ji _i = n, T ; 

T 2 ,i + 1, otherwise 



2,i - 



= n-l 



and, 




1 




1 



finally, 
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Then i k k +l {T) = V. Conversely, let V be the inverse image of T e B((k + 1)A 2 ), 
if it exists. If so, and if V e B(kA 2 ), then i k k+l {T) = V, otherwise <| +1 (T) = 0. 
Finally, 

f^o-o^o^T), ifl<J<k<S 

i%T) = I i k k+1 o . . . o t j : ; o L y\T), if 1 < k < j < s 

[T, if 1 < j = k < s. 

The partially ordered set of branching components of B(kA 2 ) also has a sym- 
metry with respect to the rank function relative to s, called *£>C-duality. If B is 
a branching component of B(kA 2 ) with rk(B) = r, then there is an isomorphic 
branching component B* BC with rk(B* BC ) = 2s — r which is the image as a set of B 
under the Lusztig automorphism [23(. For T 6 B C B(kA 2 ), we define T* BC to be 
the element corresponding to T in B* BC . Then *BC commutes with L k for all j for 
which L k (T) 0. Using this fact, we may compute T* BC explicitly as follows. Let 
I = min{j|4(T) ^ 0}. We define a map ip in B(IA 2 ) as follows. For T 6 B(IA), 
let a = ^{i\T\,i = l},b = #{i|T 2j j = 1}. Iteratively compute the following: for 
1 < i < a — 1 let T' Xi = T\^T' 2i = T 2ii , and let T 2a = T 2ja . Starting from i = a 
until we reach i = m such that Ti, m +i > T 2m and (Ti >m+ i, T 2 m ) ^ (n,n), (n,n), 
and (T l m+1 , T 2i?n+1 , T 2 m ) 7^ (x,x,y) with x < y < a; , if such an index, m, exists, 
or I — b otherwise: 



if ^~2,i+l 7^ 

T u = ^ " - 1> if T 2,i = n, T M+1 = T 2ji+1 = n 
otherwise, 




if T 2 ,i+i 7^ 



M+l 



T 2,i+i = { n - 1, if T' 2i = n, T M+ i = T 2)t+1 = n } ,i^m 

otherwise, 



and let 



T 2t m, if T 2jTn 7^ Tj 



x = { n - 1, if T 2jm _i = n, Ti, m = T 2>m = n 
^2™ — 1, otherwise. 



Then, let 



if x ^ T 2>m+ i, or m = I 
T i, m = { n, if T ljm+1 = n, T 2 , m+1 = x = n - 1 
x + 1, otherwise, 



and, 



,m+l ; if x ^ T 2 

,m+l 

T 2, m = { n, if T hm+1 = n, T 2 , m+1 =x=n-l } ,m^l 

T 2 , m +i + 1, otherwise, 
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and, for m + l<i<l — b — 1 



T' 

1 2,i 



n, if T M+ i = n, T 2ji+1 

Ti i + 1, otherwise, 

n. 



n 



if 7^ ^2,i+l 

if T M+1 = n, T 2)i+ i = 4 = n 



T 2ii+ i + 1, otherwise, 



Finally, let T' 2 l _ h = 1, and • = ly in all other cases. Defining t/>(T) = T", we 
have: T* BC = \^ a - b {i l k {T)). 

Finally, for T e B(k A2) C -B 2,s we have the following: 



a(T) 


_ L s+l-k^rp*BC 


e T 






= °CfMT)) 


£o(T) 


= e,{a{T)) 







where / = min{j|4(T* BC ) ^ 0}. 
Example: Let T = 



be an element of the £>f ^crystal B 2 ' 2 . T* BC =. We 



compute eo(T). The minimum I G {0, 1, 2} such that 4 CO 7^ is 1, with 4 CO 

and we have tr(T) = i?T* BC = 4 1 



1 




1 




2 


. Then r/> ^ 




- 




2 


2 




T 



( 


2 




2 




T 


- 


T 



Then 



On this tableau, we have (ei(er(T))) 



*BC 



and the minimum 



V e {0, 1,2} such that 4'((ei((r(T)))* BC ) ^ is 0, so we have e (T) = (x(eier(T)) 
2 



D. Lecouvey's insertion algorithm for D n -crystals 



Let 61 ® 6 2 <8> &3 ® • ■ ■ <8> 6; be a highest weight element of the D n -crystal B(Ai)® 1 
and B{pi b 2 ... 6/) be the _D n -crystal generated by bi <g> b 2 <S> &3 <S> ■ ■ ■ ® h. We 
review the "column insertion" algorithm in 21] . Let £ : 5(1 2 1) 



5(1 1 2), be 
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the unique isomorphism between the given crystals. Explicitly: 



x 
x 



V, 
z, 



if z < x < y,y ^ z, 
if x < z < y, y ^ x, 
(z + 1) <g z + 1, if y = z, z < n — 1, z < x < z, 

xA<x<n,x<z<x, 



1 ( 
n 
n 



1) <g> z, if (x,y,z) 
if (x,y,z) 
n, if (x,y,z) 



x — 1 <g) (x — 1) <g) z, if y 

y <S> x <S> z, if y > n — 1, (x, z) — (n, n) or (n, n) 

x ® z <S>y, if z < n — 1 , (x, y) = (n, n) or (n, n) 

n, 



Let 



-l 



[n, n, nj or (^n, n, n) 
(n, n — 1, n — 1), 
(n, n — 1, n — 1). 

5 5(1 1 2) ® ^(l)^- 2 



& : S(l)®*- 1 g> 5(1 2 1)® ^(l)^- 2 5(1) 
be the following crystal isomorphism: 

We use these maps to define the insertion of a letter b G B into a column T, notated 
6 — > T, of a Kashiwara-Nakashima tableau ([l8j]). Let T = Tii <g> T 2j i <g • • • <g) Tj^i, 
where Tn E B,l < i < I be the reading of T. 



1. If T = then 6 ->• T = [&]. 



2. ifT 

6^ 



Tn 5 m other words if it consists of a single box and 6 < Tn, then 



Tl,l 




6 


^1,1 



3. if b > Tf- t i or (Tj(. i, 6) = («-,?!), or (n, n), and there exists some y < n G i3 in 
the sequence S = (T lil; T 2) i, . . . , T fc l , 6) such that y is also in S and #{x G 
S^x < ?/ or x > >?/, then b — > T is the column formed by removing (z, 1) 
where z is the least such letter that occurs, 

4. otherwise, if b > T fc l or (T fc l ,6) = (n,n), or (n, n) then 6 — >• T is just the 
column T with [foj appended to the bottom, 

5. otherwise, it is the case that b < Ti^, k > 2. In this case, let T[ 2 ®T[ X ® 
TLi---®T[ k = 66 • • • &-i(Ti,i <g> T 2 ,i ® • • • ® (8) 6). Then: 



T' 


T" 


T' 

i 2,l 






7" 
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Now, let T = T X T 2 ■ ■ ■ T k , where T 1 , T 2 , . . . , T k are the columns of T . Define the 
insertion b — > T, b G B recursively as follows: 

1. Compute b — > T 1 as above. If case 4 occurs, then 



2. If case 3 occurs, let b — > T = T[ x ® T[ 2 <8> • • • ® T[ t for some / > 0. Then 
b — > T is given by successively inserting the letters of (T 1 )' into T\T l as 
follows: 



3. If case 2 or 5 occurs, then notice that in both cases, b — > T 1 = (T 1 )'(T 2 )' 



where (T 1 )' is a column and (T 2 )' = c , for some c E B. Then we define: 



Lecouvey |21( defines an oscillating tableau Q of type D n to be a sequence 
Qo> Qii Q2, ■ ■ ■ 1 Qi of pairs Qk = (Ok,£k), where Ok is a Young diagram whose 
columns have height < n and G {— , 0, +}, satisfying, for k — 1, 2, . . . , I: 

1. Ok and Ok+i differ by exactly one box, 

2. e k ^ and e k+1 ^ imply e k = e k+1 , and 

3. e k = if and only if O k has no columns of height n. 

Then, we have the following analogue of the Robinson-Schensted correspondence 
for type D n : 

Theorem 3 ((2l[). There is a bijection between the set B(A\)® 1 and the set of all 
pairs (P, Q) where Q = (Q = (0, 0), Qi, Q 2 , ■ ■ ■ , Qi) is an oscillating tableau, and 
P is a Kashiwara-Nakashima tableau of shape 0\ such that, if P has a column of 
height n, whose kth entry is n (resp. n) then n — k is even (resp. odd) if Ei = +, 
and n — k odd (resp. even) if E\ = —. 

Explicitly, this bijection is given by sending 6i®&2®- • -®h G B(ki)® 1 to (Pi, Q), 
where P k = b k ->• -+ (• • - 61 ->• 0)), 1 < k < I, Q k = (O k , e k ), 1 < k < I, O k 
the shape of the tableau P k , and e k = is P k has no columns of height n, e k = + 
if P k has a column of height n whose /cth entry is n (resp. n) such that n — k is 
even (resp. odd), and e k = — otherwise. 

Example: The element 2®4cg>4cg>3of the /^-crystal B(Ai)® A corresponds to 
the following sequence of tableaux: 



b^T = (b^T l )T 2 ■■■T k . 



b^T = T[ A (T/_ X)1 -> ( > T[ x T 2 T 3 ■ ■ ■ T k )). 



b^T= (T 1 )'(b^T 2 T 3 ---T k ). 



12 



2 -> = g 

*->@=lf. 



4^ 



_2_ 

7 

4 



3 -> 



. We compute £i£ 2 (2®4®4® 3) = £i(2®4®3®4) = 4 ® 2 g> 3 ® 4 



^12 ® ^i i ® ^2 1 ® ^3 1- The resulting tableau is: 



E. Combinatorial R : B r > s ® £ r '- s ' -> £ r '' s ' ® £ r '< s 

In this section, we define the combinatorial R matrix for B r,s ® B r ' ,s ' and the 
energy function H : B r ' s ® £T ,s — > Z. 

Proposition 1 ([15]). There exist a unique crystal isomorphism ft : B r,s ®B r ' ,s ' — >■ 
5 r > s ®£ r - s and a function H : B r > s ®B r > s -> Z unique up to an additive constant 
satisfying the following property: for any b G 5 r ' s ; 6' e5 r - s , andi e I such that 
ei{b®U) ^ 0, 

r ff(6 ® 6') + 1 t/i = 0, <^o(fe) > £o(6'), ¥>o(6') > e„(&), 
H{ei{h®b'))= I H(b®b')-l if i = 0, ^{b) < e (V),<p (V) < e (6), (10) 
^ i7(6 <g> 6') otherwise, 

where V ®b = lZ(b ® &'). is called an energy function on 5 r,s ® B r '' s ' . 

Proposition 2 (Yang-Baxter Equation). The following equation holds on B r,s £g> 
5 r '' s ' ® B r "' s " 



(ft ® 1)(1 ® ft) (ft <g> 1) = (1 ® ft) (ft ® 1)(1 <8> ft), 
where 1 denotes the identity map. 
We define ft^+i to be the map 



'11' 



I®* -1 g> ft <g> l®^- 1 : (g) ->■ (g) E^M-Mfc) 

fc=i fc=i 
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where r transposes i, i + 1 and fixes all other integers. In this notation, (ITTj) 
becomes: 

(12) 

for L > 3, z = 1, 2, . . . , L — 2. 

Definition 4 (Affinization). T7ie affinization of B r,s is defined to be the set 

{z n b\n eZ,be B r ' s } 

where the action of e,, /j is defined as: 

Bi{z n b) = z n+5M ei{b) 
f t (z n b) = (6). 

The combinatorial iZ-matrix ft Aff : Aff(£ r ' s ) <g> Aff( J B r '' s ') -»■ Aff( J B r '' s ') <g> 
Aff(£>' r ' s ) is given by: 

where 6' <g) 6 = K(b ® b'). 



F. Combinatorial .R-matrix for B 1,s ® B 1 ' 8 ' 

The combinatorial .R-matrix for the .D^-crystals B 1,s £g> 5 1,s ' has been given in 
[?J • For the tableaux: 

b®b' = (x 1 ,x 2 ,...,x 1 )(g)(y 1 ,y2,...,y 1 ) E B 1 ' 8 ® B 1 ' 8 ' 
we set z = min(xi, y±) so that b®b' are associated with the tableaux: 



1 


1 




1 


T 





V2 




Vk 


1 


1 




1 





where k = s' — z. Now, define the tableau — v% — >■ (v 2 —>■(-•• Vk — >■ T*) 
which is necessarily of the form, for I = s — z\ 



Jl 




3k 






k 


h 




I'm 





for some < m < L^J- Next, we use reverse column insertion to remove the 



boxes containing ij, 



!)•••) f m+l j 'mi j 



zi in order. Let Wi, w 2 , ■ ■ ■ , wi be the 
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sequence of letters that are produced in each step, and T^\T^ 2 \ . . .T^ denote 
the sequence of tableaux. Then 



K(b ® b') 



1 


1 




1 


y(0 





w 2 




Wl 


1 


1 




1 







The energy function is given by: 

H (b ®b') = 2 min(7, k)-m- 2s 2 - 



G. Combinatorial i?-matrix for A n 

Define B r > s to be the set of r x s semistandard tableaux in the alphabet 
jl',2', . . . ,n'}. The crystal structure of this set is given in j29[, and the combi- 
natorial R matrix is given as well in loc. cit. (see also 34|) by the following row 
insertion procedure. 

To an element 



T 



Ti,i 


Tl,2 




Ti,s 


T2,l 


T2.2 




T2,s 












T r ,2 




T 

- 1 r,s 



e B 



r, s 



we associate the row word row(T) = T\ iS (g)Ti iS _x <S>Ti^ <g> T 2iS (g)T 2iS _i <S> ■ ■ ■ ®T 2t i ® 
■ ■ ■ <g> T r ^ s ® T r)S _i ® ■ ■ ■ ® T n i. Then we have the following, which follows from the 
fact that the decomposition of B r,s ® B r ' ,s> is outer-multiplicity free as v4 n -crystals: 

Theorem 4 (@, see Q). Let T ®V E B r ' s g> B r '' s ' . Then: 

1. the map U : B r ' s ® B r ' s -> B r ' s ® B r ' s is given by the following: 

K(T ® T') = V ® f 

if and only if row(T r ) — > T = row(T) — > T', where row(T) — > T' denotes the 
resulting tableau from row inserting all the letters of row(T) in order into T' , 

2. let d(T, T') denote the number of nodes in the shape of row(T) — > T' that 
are strictly to the right of the max(s, s') column. Then the energy function 
H : B r ' s ® B r '' s ' — > Z is given by 



H(T ® T') = d(T, T') — min(r, r') min(s, s'). 
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To the tableau T e B r,s we associate the coordinatization 



f Xl^l X\2 ■ ■ ■ 3^1,n\ 
3^2,1 x 2,2 ' ' ' ^2,n 



\X 

r.l %r,2 



fc ^ >Q . 



III. COMBINATORIAL i?- MATRIX FOR D 



(i) 



In this section, we give the combinatorial i?-matrix for the U q (Dn )-crystals 
B 2s ® 5 2,1 , first for the highest weight vectors and then for arbitrary vectors using 
the analogue of the Robinson-Shensted correspondence (Theorem [3]) • 

The D n , n > 4 highest weight vectors of B 2 ' s ® B 2 ' 1 £1X6 cLS follows: 



Highest weight vector 



u k ® 



< k < s 

u k <S> 0, < k < s 

1 < k < s 



U k <8> 



,1 < < S 
, 1 < k < S 

,2 < k < s 
, 1 < k < s 
, 1 < k < s 
,n = 4,1 < k < s 



Classical Weight 



(A; + 1)A 2 
kA 2 

( Ai + {k- 1)A 2 + A 3 ,n > 4 
\ Ai + (A;-1)A2 + A3 + A4,n 
(k - 1)A 2 + A 4 ,n > 5 
(k - 1)A 2 + A 4 + A 5 ,n = 5 
(k - 1)A 2 + 2A 4 , n = 4 

2Ai + (fc - 1)A 2 



Ai + (k- 2)A 2 + A 3 ,n > 4 

Ai + (k - 2)A 2 + A 3 + A 4 , n = 4 



A;A 2 

(k - 1)A 2 

(A; - 1)A 2 + 2A 3 



We will also use the D n , n > 4 highest weight vectors of B 1 ' 1 ®B 2,1 and B 2,1 ^>B 1 ' 1 : 
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Highest weight vector in 
B 1 ' 1 <g> B 2 ' 1 



a 











Classical 
weight 



Ai + A 2 

A 3 + <5„ j4 A 4 

Ai 

Ai 



Highest weight vector in 







5 




Classical 
weight 



Ax + A 2 

A 3 + <5„ i4 A 4 

A x 

A x 



A. Combinatorial i? matrix for highest weight elements 

We begin by computing the 0-string through certain elements of B 2 ' s . 



Lemma 1. If T 



T2,l 



given by the following: 



G -B(A 2 ) C B > s , s > 1, then the 0-string through T is 



fF = { 



v- 


-1 


1 




-1 


Ti 


,i 


v- 


-1 


2 




-1 




v- 


-1 








-1 


1 




V 


T 


i.i 




1? 


T2,l 


i 



o, 



tfT 2tl =%T^i{l,2},l< 3 <s 
ifT 2:1 =T,T hl £{2,2},l<j<s 
ifT-- 



,!<3<s + l 



(14) 



otherwise, if 1 < j < s — 1 
otherwise. 
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Similarly, we have: 



eiT = < 





T 


-1 


I 


V 


- 1 


T2,l 


T 


-i 


2 


v 


-i 








v- 1 


•> 




Ti,i 


2 3 




T2,l 


V 





i/T 1)1 = 2,T 2)1 ^2,l<j<s 
i/T 1)1 = l,T 2)1 ^2,l<j<s 



1 < j < s + 1 



(15) 



0, 



Also, the 0-string through T 



1 


2 


3 


T 



i/T 2)1 £{l,2},l< j<s-l 
otherwise. 

is given by: 



? T 



p-l 


1 


3 


2J'-i 


3 


3 



1 




3 





, 1 < J < s - 1 
1 < j < s - 2 



Proo/. For T 



T 2 ,i 



e B{A 2 ) C B 2 > s , s > 1, we compute /qT = (jo/ 1 o ^'(T) 



a(f(a(T)), where we have used the fact that cr = <r 1 . In all cases except T = 

' ^(A 2 )): 



For all T G I?(A2), the following may easily be verified (recall that 



rpf rp 



*BC 



T 2 a 



, ifT 2jl = l, 



^ T, otherwise. 
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Then, we have: 



a(T) = l\ +1 - 1 {T') = { 



K-i 



T 2 a 



M-i, 



if T 1A = 1 
ifT 2il =T 
if T - 



T 2 a 



s-2 



T 



, otherwise, if s > 1 
otherwise. 



where j\fh 



1 L*/2J 


2^ (mod 2) 


2 Lfc/2j 




2^' (mod 2) 


T Lfc/2j 



denotes the null configuration of size k > 



[28|). We compute: 

K-i-j 



7? 


T 2 ,i 


V 


1 



Ti, 



Ti,i 



2-> 



fi a (T) = 



K-l-i 



2-' 



if Ti.1 = 1, 1 < j < s - 1 
ifT 2jl =T,l<j< S -l 
ifT 2jl =T,T ljl ^2,j = S 
ifT 2A = T,T 1A = 2,j = s + l 

l < j < s - 1 



if T 



1 




2* 


Ti.i 


22,1 




I 


2 


2 s-2 


^1,1 


5 


^2,1 


yS — 2 


1 



2«-i 

-rs-l 



Ti.i 



otherwise, if s > 1, 1 < j < s — 2 

otherwise, if s > 1, 1 < j — s — 1 

otherwise, if s > 1, T 2j i = 2, 1 < j = s 

otherwise, if s = j = 1, T 21 = 2 
otherwise. 



19 



Thus we have: 



,i'-i+i fj 





T2,l 


V 


1 



2-> 



i-l 



r 1 



2J'-i 



2- ? 



2^ 



UT hl = l,l< 3 <s-l 
if T 2)1 = 1,71,! ^2,1 < j <s 
if T 2il = 1,71,! = 2,1 < j < s + 1 
if 71,! = 2,72,1 =T,1< j<s-l 
■ , 1 < J < s - 1 



if T 



2 


2i-2 




72,1 




T 


2 s-l 






T- 1 


i 



u 



otherwise, if s > 1, 1 < j < s — 1 

otherwise, if s > 1, T 2 ,i = 2, 1 < j = s 

otherwise, if s — j — 1, T 2 ,i = 2 
otherwise, 



where /' = rcim.{j\L 3 s _ 1+ J{a(T) ^ 0}. Applying *i3C on each of the tableaux gives 
the desired result. A similar computation gives the result for e Q . 



For the 0-string through T 



1 


2 


3 


T 



we have I = 2, and compute: 
1 



fii s 2 (T* BC ) = { 



1 






2? 


2 


y 


-2 


V 


I 


2 


2 s-2 


2 






3 


V -2 


T 


1 





2 


2J'-i 


2 


3 


P- 1 


T 



o, 



if 1 < j < s - 2 

j = s - 1 
otherwise 

1 < j < s - 1 
otherwise 
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(since I' — j + 1 in all cases) 

(i i :(R4(t* bc W bc = 



p-1 


1 


3 


2 j-i 


3 


3 



o, 



1 < j < s - 1 
otherwise 



A similar computation gives the result for e 3 . 



□ 



Now we are ready to prove the following: 
Theorem 5. On the D n highest weight vectors in B 2 ' s ® B 2 ' 1 , s > 1, we have: 



K(u k <g> b) = < 



<8> Mfc, 



Mi 



Mi <g) M fe _ 2 , 



Ml 





3 


2 fc-l 


4 



ifb 
ifb 
ifb 



Mi (g) 


Ms, 




if b — 


1 1 -i 


® M s , 




ifb = 


Ui 


Mi <g) 




2 




ifb = 


Ui 


2 fc+l 


I 


• 


Mi (g) 


p-1 


1 




ifb = 


1 


2 s-l 


3 




3 


Ml ® 




1 


3 


, ifb = 


1 


2 fc-l 


3 


3 


3 


Mi (g) 




3 




ifb = 


3 


2 k-l 


4 




4 


Mi <g 




3 




ifb = 


3 


2 k-l 


3 




3 


Mi <g 




1 




ifb = 


1 


2 k-l 


2 


• 


2 


Mi (g 


^k-2 


1 




ifb = 


3 


2 k-2 


3 




2 



and k = s, 
and 1 < k < s — 1, 
and 1 < k < s, 
and 1 < k < s, 
and 1 < k < s, 
and 2 < k < s, 



if b = and k = s, 

ifb = andO <k < s - 1, 



and = 1, 

and 2 < /c < s, 

and 1 < k < s,n = 4. 



(16) 
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Proof. Since 1Z is a Z?n crystal isomorphism, we have wt(6 <8> b') = wt(lZ(b <g> 6')). 
Since u s ®Ui is the unique highest weight vector of B 2,S ®B 2,1 such that wt(&<8>&') = 
(s + 1)A 2 , we must have 7£(w s <8> u%) — u% <g> w s . Therefore lZ{e 3 {u s <8> Mi)) = 
e 3 TZ(u s <8> Mi) = 6q(mi Cg> w s ), j > 0. By the previous Lemma, and Definition [2J 



m s <S> 0, 



Us-j+2 



if J = 1 

if 2 < j < s + 2 

if s + 3 < j < 2s + 2 



and, 



g) U a ) 



Ml ® U s -j, if 1 < j ' < s 

, if s + 1 < j < 2s 



0® 



T 


-s 


T 


-s 


2 s 




V 


) 


2 s 




I s 


• 



if j = 2s + 1 
if j = 2s + 2 



Setting k — s — j + 1, we see: 

lZ{u s ® 0) — Ui <g> M fl _i, 
and, setting k = s — j + 2 gives: 



7?. M fc 



«i ® Ufe_2, if 2 < < s 
U\ <g> \^=\ , if = 1 



Also, we have 
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however, these are not D n highest weight vectors. By using column insertion 
(Theorem [3D we find the corresponding D n highest weight vectors, and obtain: 



TZ{u k ®u x ) = Ui ® 





2 


2 fc+l 


I 



, < k < s — 2 (and nothing else new). 



Ti, 



T2A 



or 



4 



Now, we consider the case £g> T, where T 
if n = 4. We see that the corresponding highest weight vector is <g> u\. We have 
seen that 1Z(0 g) Mi 



Ml ® 



the tableau T): 7^.(0 (giT) = U\ 



which gives (using inverse column insertion on 
Acting on both sides by Jo a sufficient 





2 


^2,1 


T 



number of times using the previous Lemma, and commuting with 1Z gives: 

if T - 



n(u k ® r) 



Ml 



Ml 



Mi ® 



Mi (g) 



Ml 



Mi 



Mi (g> 



s-1 



)S-1 



fc-1 



)fc-l 



i fc-1 



)fc-l 



I fc-1 



)fc-l 



fc-1 



)fc-l 



I fc-2 



)fc-2 



I fc-1 



)fc-l 



if T 
if T 
if T 
if T 
if T 
if T 



and k = s, 
and 1 < k < s — 1, 
and 1 < k < s, 
and 1 < k < s, 
and 1 < k < s, 
and 2 < k < s, 
and 1 < k < s, n = 4. 



Finally, consider 



0. The corresponding highest weight vector is u s <8> 0- We 







2 




2 s_i 


T 




r- 1 



have seen that lZ(u s <S> 0) — u\ <S> m s _i, which gives 1Z 
Acting on both sides by fo sufficiently many times, and commuting with TZ gives: 

K(u k ®T) = 



<g> M fc , if T = 0, < < s - 1 
<E> m s , if T = Mi, fc = s — 1. 



□ 
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From the computation in the proof of Theorem one may easily deduce the 
following: 



Corollary 1. For the highest weight vector Uk <8> T e B 2 ' s <S> B 2 ' 1 , < k < s, 



H(u k ® T) = < 



0, ifT 
-1, ifT 
T = 





1 




2 




1 




2 


1 


,k 




3 





k = s 



k = s — 1, 



,k = s, or, 



T = } k = s, 
—2, otherwise. 



(17) 



Similarly, we have: 



Theorem 6. On the D n highest weight vectors in B > ® B ' , 



n(^\®T 



0, ifT 
0, »/T 







0, ^/T 



Proof. The cases where T 



1 




2 






or 




are 


2 




3 





are immediate because there is no outer 
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multiplicity in these cases. Consider 7?- T [T] <g) 0) . We have: 

fo (0 «)0) =0 ® 

; ^(g ®0) = | ® g 

ft([T] (g)0) =e I ^ <8) 



Finally, the case where T 
possibilities. 



is immediate, since we have ruled out all other 

□ 



B. Combinatorial R matrix for B 2,s (g> B 2,1 

Recall (Theorem [3]) that there is a bijection between the set B(Ai)® 1 and the 
set of all pairs (P,Q) where Q = (Qq = (0, 0), Qi, Q 2 , ■ ■ ■ , Qi) is an oscillating 
tableau, and P is a Kashiwara-Nakashima tableau of shape Oi. In fact ( [21] ) , the 
oscillating tableau Q enables us to determine the highest- weight vector in B(Ai)® 1 
corresponding to the pair (P, Q): namely it is the vector b\ <S> b 2 <8> • • • <S> bi, where 



n. 



n. 



if 1 < % < n, Ok has one more box in the ith row than Ok-i 
if 1 < i < n, Ok has one fewer box in the ith row than Ok-\ 
if Ek = +, 0(resp. — ) and Ok has 1 more (resp. fewer) box 
in the nth row than Ok-i 
otherwise. 



The following procedure then allows us to compute the combinatorial i?-matrix 
for any element T ®T' e B 2 ' s ® B 2 ' 1 . 

1. Insert the letters of T' from top to bottom into T, keeping track of the rows 
in which boxes are being added or removed and the positions of any n or n 
in a column of height n. If there T' = then do nothing. Call the resulting 
tableau P. 



2. Find the highest weight vector b\ ® b 2 <8> ■ ■ • from Step 1. The result gives 
the reading of a unique highest weight vector Uk <S> T" G B 2,s ® B 2 ' 1 . 
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3. Compute TZ(uk (g T"), and interpret the result as an element of £?(Ai)® J , 
and use it to deterine the oscillating tableau Q'- 

4. For the pair (P, Q'), reverse Lecouvey's algorithm to get a sequence of letters. 
The resulting sequence of letters gives the reading of a unique pair of tableaux 
T ® f G B 2 ' 1 g> B 2 > s which is 1Z(T g> T"). 



Example: We compute 7?. 



as follows: T" 



2 1 



4 




1 


4 


<g 


2 


4 


)- 


4 





for the -D 4 -crystal £? 2,2 <g -B 2 ' 1 . We insert 



1 


4 


2 


4 



. The corresponding highest weight 



vector is 1 <g 2 <g 1 <g 2, which corresponds to Mi ®«i G 5 2 ' 2 <g -B 2 ' 1 . By Theorem 
[5] we have 7^(iti <g u\) 



1 


1 


2 


2 



M2 = <g 

tableau T' successively from rows 2, 1, 2, 1 which gives 4 <g 4 <g 1 <g 2, which, in the 
B(0) <g £?(2A 2 ) component of 5 2 ' 1 (g B 2,2 is interpreted as 



. So, we remove letters from the 



1 


4 


2 


4 



IV. SOLITON CELLULAR AUTOMATON 

We define V L = {h®b 2 ® ■ ■ • ®b L e (B 2 > l )® L \T n = m, for n sufficiently large} 
to be the set of states ofthe Dn soliton cellular automaton. We depict the 
operation TZ(b (g b') = b (g b by: 



b'- 



Fix s > and let u s G B 2,s be the highest weight vector for D n . For p = 
b\ (g b 2 <E> ■ • • <g &l G Pl we define the izme evolution operator Ti(p): 

r,(p) ® U , = n LL+1 n L ^ L ■ ■ ■ n 23 n 12 ( Ul <g> P ). (19) 

The transition of phase <g 6 2 <g • • • <g = bi ® b 2 ® ■ ■ ■ ® b L is depicted as 



. M (2)... M (L-1) 



&i 
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We define the state energy to be the sum: 

L-l 



(20) 



i=0 



We can also use the combinatorial .R-matrix TZ : B 1,l ®B 2 ' 1 — > B 2,1 ®B 1,1 to define 
an operator Th similar to TJ, by 



T\ (p) ® 6(p) = l ' ' ' n 23 n 12 ( [T] <8) 



(21) 



In this case, 6(p) is dependent on the state p so we indicate this dependence in the 
definition. 



A. D^-solitons and their scattering rules 

Experience from other soliton cellular automata has shown that states p G Vl 
satisfying Ei{p) = 1 correspond to the so-called "one-soliton states". 



Proposition 3. In the Dn^ SCA, Ei(p) = 1 if and only if p ^ uf L and is of the 
following form: 

wf*® 



2 




2 






02 



2 




1 




1 


bj 




bj+i 




bj+2 



® uf 



(22) 



for some i,j<k,l& Z> such that i + k + / = L, and some bi > b 2 > • • ■ > b k G 
B\{1,2,2,I}. 

Proof. Let p = 7\ ® T 2 ® • • ■ g) T L G P L be such that £ a (p) = 1. For T ® V G 
<g> 5 2 ' 1 it is the case that K(T ® T') = T ® T'. Thus, £?i(p) = <8> 
Ti) — X^=i -^(^i ® ^i+i) — 1) an d & H t ne terms appearing in this sum are > 

by Corollary [TJ Now, u\ <E> T = if and only if T = 



. Suppose that we 



have H(T k ® T fc+1 ) = —1 for some T k = U\. In our SCA, T k = Ui,k S> and 
if(0 g) T') < 0,T' G -B 2 ' 1 , hence it must not be the case that T k+1 = 0. So 



we are left in the case that u\®T has D n highest weight vector u\ 

column insertion, (Theorem [3]) this is only the case if T 

b G B\{1,2,2,I}. 

Since £d(p) = 1, we must have <8> 



1 






, or 


6 





By 
where 



be of the form 



bx 



0,k < j < L. Supposing Tj to 
b\ G {1,2}, we see that T J+ i must be 



b[ 



b' 



,b[ < h,b' 2 < b 2 ,b' 2 G 
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B\{1, 2, 1} in order for Tj®Tj + \ to have D n highest weight vector u\ <g) u\. Finally. 
- ,j > for our SCA. 



□ 



The proof of the following proposition is completely analogous to that of Propo- 
sition 9 in j26[, except that £i(itj), and (pi(u{) are now 0. 

Proposition 4. Le£ p G "Pl- If e~i(p) ^ i/ien eiTi(p) = 7}(ej(p)),z 7^ 0,2 and 
Ei(ei(p)) = Ei{p), otherwise ejTj(p) = 0,2 7^ 0,2. Tae same relations hold for 
/i,«^0,2. 

Recall the -crystal _B 1,S = {(xi,x 2 ) G Z> |a;i + x 2 = s}, where (xi,x 2 ) can 
be associated with the set of tableaux: 



V 


V 




V 


2' 


2' 




2' 








X\ 






X2 





and, similarly, the Dn ''-crystal B 1,s = {(x\, x 2 , • • • , ^1) G 1%?q\X\ + x 2 + ■ ■ ■ + x\ 
s, x n = or x n = 0} can be associated with the tableaux: 



1 


1 




1 


2 


2 




2 




1 


1 




1 
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Proposition 5. Let (T, T") 6e an element of the ®D n Z 2 - crystal B lyS x i? 1,s , n > 
5. Define the map i s : B l ' s x £> 1,s — >■ (i? 2 - 1 )® 5 owen fry 



( l'i 2 /s_J ' 



61 


&2 




6 S 



2 




2 









b[ 



61 



61 



-3-2 



61 



where, 



b' 



b k + 2, ifl<b k <n-2 



c + 2, if bk = c, I < c < n — 2 
Then i s satisfies the relations: 

i s (e?T,T') = e 1 i,(T,T f ),i a {ffT,T r ) = A^(T,T'), 

i s (T,~e?T') = e i+2is (T,T'),z s (T,f~ t D T') = f i+2 i s (T,T'), l<^<n-2, 
If n = 5 taen we have: 



3' 



2' 



1' 








2' 








3' 








1 — ^ 


5 






1 — y 


4 




4' 


1 — y 


3 


4' 






1 


4' 






5 
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yields a map i s satisfying 

i s (e?T,T') = e x is{T,T'),i s {fiT,T') = fxi s (T,T'), 

i,(T,e?r) = e c{i)+2 i s {T,T'),i s {Tj*T') = £ (0+2 i fl (T, T'), 1 < i < 3 

where a transposes 1 and 2 and fixes 3. 
For n = 4 we have: 

((xi, x 2 ), (j/i, y 2 )) 3 s_max(a:i ' 2/l) (8) 4 ( " l)+ ® 4 ( ^' l)+ ® 3 min(xi ' yi) 

where (x) + = max(x, 0), yields a map i s satisfying: 

i s {e?T,T',T") = eMT,T',T"),z s (ffT,T',T") = f\i s (T,T' ,T"), 
i s (T,~efT',T") = e 3 i s {T,T',T"),i s (T,ffT',T") = f 3 i s (T,T' ,T"), 
i s (T, T', efT") = e A i s (T, T', T"),i s (T, T', f*T") = / 4 z a (T, T', T"). 

Proof. We have: 

,T') = i s (e?(2® s - 3 



and 



e x i s ( l' j 2' s ~ j ,T') 



ei 



2 




2 


6i 




Ci 



2 




1 




1 


6' ■ 

s — 1 




fc' -1 







2 /s ^' ,T'), si 



since 



Observe ei only has non-zero action on the top row of i s ( l'i 
b' k G {3, 4, . . . , n, n, . . . , 4, 3}. Also, the action of e\ on the top row of i s ( ].'•? 2' s ~i ) 
is equivalent to ei{2'® s ~^ ® 1'® J '), from which we may deduce that i s {e^T,T') = 
eii s (T,T'). The proof is similar for i s (f^T,T') = fii s (T,T'). 

The proof of i s (T, ~e?T) = ~e, l+2 i s {T, T% and i,(T, f\ D T') = f l+2 i s (T, T"), 1 < 
i < n — 2 is similar, except that , e~i +2 and fi +2 act on the bottom row of i s (T, T'). 
One easily checks that the actions of e~i +2 and fi +2 on b' k are equivalent to ef and 
j?on& fe . □ 

Remark: The above map is a bijection from I? 1 ' 5 x B 1,s , (or i? 1 '* x B 2,s , or 
(fi 1 ' 3 ) 3 ) to {p G V L \Ei(p) = 1} by Proposition! 

A state of the following form is called and m-soliton state: 

• • • [si] [s 2 ] [s m ] . . . 

where s\ > s 2 > ■ ■ ■ > s m , . . . [s] . . . denotes a local configuration of the form 

for some b\ > b 2 > 



2 




2 


bx 




b 2 



® •••® 



B\{1,2,2,1}, and the [s* 



1 



i 



6; 



'i+2 



1 



> fes G 



are separated by sufficiently many «i's. 
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Proposition 6. Let p be a one-soliton state of length s. Then 

1. E k {p) = min(k, s), 

2. T k {p) is obtained by rightward shift by mm(k, s) lattice steps. 



Proof. By applying sufficient operators e iy i ^ 0,2, p becomes 
■ ■ ■ ®U\. Using the combinatorial R- matrix we obtain: 



K 



V 


■yl— i 


2* 


3 s-i 



Ui 



K 



K 



V 




2* 


3 s-i 



v- 


i 


■ys—i+l 


T- 


i 




I s 






3 s 






p+i 


^ s — i — 1 


2 i+i 


y-i-1 



if % > 0, 



if i < s, 



1Z(u s <S> Ui) = Ui <S> u s 



and 



H 



r 


-^s — i 


T 


3 s-i 



H 



H 



V 


ys—i 


T 





<8> Ui 



-1 if i > 0, 

0, 
0. 



Ui ® Ui 



(23) 

(24) 

(25) 
(26) 

(27) 
(28) 
(29) 



where, as before, the symbol V means that b is repeated j times. If k < s then 





1 


(^Uk <S> 






3 



<g> uf L - s 



= Ui 



and 



L-l 



s fc (p) = -X)^( u(<) ® 6 <) = fc - 



i=0 
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Otherwise, 



T k (p) ®U s =K LL+ i--- 7^23^-12 (Uk ® 



<8> uf L - s ) 



u 



®L-2s 
1 



® Uk, 



and 



L-l 



i=0 

The result follows because the z 7^ 0, 2 commute with T s and preserve £7 S (Propo- 
sition @J. □ 

We now consider the two-soliton case 

p= ... [ Sl ] . . . [s 2 ] . . . 

where si > s 2 . We can use Proposition[3]to associate a two-soliton state T*(p) := p t 
at time t with the element z~ kl b x ® z~ k2 b 2 e AS(B^ Sl x B 1 ' 81 ) ® AS(B l ' S2 x 5 1 ' 52 ), 
where fcj := — min(r, Sj)t + 7^, where 7, is the number of 'iti's to the left of [sj]. 
If r > S2, then we expect to see the longer soliton catch up with and collide with 
the shorter one, and, after sufficiently many time steps separate out into another 
two soliton state. At the end of this section, we will prove that is the case. In the 
following lemmas, we identify elements in Aff(i? 1 ' Sl x B 1 ' 111 ) eg) Aff(£> 1,S2 x B 1 ^ 2 ) 
with their corresponding two-soliton states. 



Lemma 2 (Analogous to Lemma 4.15 in [351 ]. Lemma 3 in |26j]). Suppose 
there is a one-soliton state p = . . . [s] . . . corresponding to z~ k ((i,s — i),b) 6 
Aff(B 1,Sl x B 1,S1 ),0 <i<s. Then T^{p) is another one-soliton state, correspond- 



ing to z k ((s,0),b), if i = s, and z k + l,s — 1 
have b{p) — if i — s and b{p) = [2] otherwise. 



l),b) otherwise. We also 



Proof. We compute Tj, 

for the two cases j = ( 

Case 1: j = 0. We compute 



2 


<8> 


2 


<g>- • -® 


2 




1 




1 


® - • -® 


1 


b s 




bs-i 


bs-j 


bs-j-i 




bs-j-2 




61 



'0 



since 67^ {1,2, 1}. Therefore, we have T^(p) = p and b(p) = [Tj . 



31 



Case 2: j > 0. We compute: 



a 0, 



K 



Jk-i 



bk-i 



where c = 1 or 2, and 6^ < 6&_i. 



□ 



Now consider the two soliton state z~ hl ((si, 0), b) £§> z~ k2 ((i, s 2 — i), b'). Thanks 
to Lemma HJ the action of has no effect on the first soliton, and therefore we 
have the following corollary. 

Corollary 2 (Analogous to Lemma 4.15 in [35[, Corollary 1 in j26|). Suppose we 
have a two-soliton state p = . . . [si] [52] • • • corresponding to z~ kl ((si, 0) x 

61) ® z~ fc2 ((i,s 2 - i) x 6 2 ) e AfftB 1 ' 31 x 5 1,Sl ) <g> AfftB 1 ' 32 x £? 1,S2 ), < i < 
s 2 . Then T^(p) is another two-soliton state, corresponding to z~ kl ((si,0),bi) <8> 
z~ k2 ((s 2 ,0),b), ifi = s 2 , and 0), 61) ® «~* a ~ 1 ((i + 1, s 2 -i - 1), 62) oi/ier- 

wzse. VFe a/so /lave 6(p) = [Tj i = Si and b{p) = \ 2j otherwise. 

Before we prove the main result on scattering of solitons, we first prove several 
Lemmas relating 1Z and T r , TL. 



Lemma 3 (Analogous to Lemma 4.17 in |35J, Lemma 4 in [26(). Assume that 
si > s 2 . For p = 2- fcl ((si,0),6!) <g> z" fc2 ((i,s 2 - i),6 2 ) e A/ft-B 1 '" 1 X B 1 '' 1 ) <g> 
AffiB 1 ' 82 x 5 1,S2 ), < i < s 2 we have 

1. T^K A ff{p))=K A B(T^p)), and 

2. b{p) = b{jl A ff{p)), 

where we use the shifted energy function H = 2s 2 + H in 1Z A ® ' . 

Proof. As in Lemma [2] we have the following two cases: i = s 2 and i < s 2 . 
Case 1: We compute 

T^K A \z- k ^ (( Sl , 0), h) <g> z- k2 ((s 2 , 0), 62))) 

= T^(z- k2+5 ((s 2 , 0), 62) ® ^""((ai, 0), 60) 

= z" fc2+5 (( S 2, 0), 62) ® ^^((si, 0), 60 



and, 



n AS (T k ((z-^(( Sl , 0), 60 ® ^ fc2 (( S2 , 0), 6 2 ))) 
= TZ AG (z- k2 ((s 2 , 0), 6 2 ) ® z-^(( Sl , 0), 60) 



z~ k2+H ((s 2 , 0), 6 2 ) ® z- k ^ H (( Sl , 0), 60 
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where 6 2 <g> 61 = 7£(6i <g> 6 2 ). 
Case 2: We compute 

T^(^ Aff (^ fcl (( Sl , 0), h) ® ^ fc2 ((z, s 2 - i), 62))) 

= T^(z- fc2+i? (( S2 , 0), 62) (g) ^""((si - s 2 - i, s 2 - i), 60) 

= z- fe2+5 ((s 2 , 0), 6 2 ) ® ^-^-^(S! - s 2 - i + 1, s 2 - z - 1), 60 

and, 

^ Aff (T h ((^- fel (( Sl) 0), 60 ® z~ k *((i, s 2 -i), 6 2 ))) 

= £ Af V fc2 ((s 2 , 0), 6 2 ) (2) ^^((z + 1, s 2 - i - 1), 60) 
= z~ k2+ii ((s 2 , 0), 6 2 ) ® ^-^^(s! - s 2 - i + 1, s 2 - i - 1), 60 
where 6 2 ® 61 = 7£(6i ® 6 2 ). □ 
The proof of the following is the same as that of Lemma 5 in [26[ . 



Lemma 4 (Analogous to Lemma 4.18 in [351 ]. Lemma 5 in [26]). Let p G Vl, I > 
0, L > 0. T6en 

1. mT l (p))=T l (T k (p)), 

2. 6(r,(p))®wi = ^(tti®6(p)). 

Now we are ready to prove the main result: 

Theorem 7. Let p = . . . [si] . . . [s 2 ] ■ ■ ■ G Vl be a two-soliton state with s± > s 2 , 
corresponding to z kl b± ® z k2 b 2 where ki,k 2 < 0. Then after sufficiently many 
applications ofT r ,r > s 2 the new state is given by 

Vl AB {z kl b x ® z k2 b 2 ) = z k ' 2 b 2 <8> z^h 

with phase shift 

k' 2 -k 2 = k 1 -k' 1 = 2l 2 + H{b x ® 62). 

Proof. By Proposition HI T r commutes with e~i,fi,i 7^ 0,2. Thus, it is enough to 
check the scattering rule for the highest weight elements 

( Aff((S 1 > s 3 ) ® Aff(( J B 1 ' S2 ) 3 ), if n = 4, 

2 fcl (( Sl ,0),6)®^ 2 ((?/!, y 2 ), b') G < Aff(S 1 ' sl x B 2 ' n )®M{B 1 ^ x B 2 ' 32 ), if n = 5, . 

[ Aff(S 1 ' sl x B 1 '" 1 ) ® Aff( J B 1 ' 52 x 5 1 ' 52 ) if n > 5. 

We will show the statement is true by induction on y 2 . 

Suppose y 2 = 0. Then j/i = s 2 . In this case, the time evolution operator is 
equivalent to that of the D^li SCA associated to the crystal B 1,s if n > 5, or the 
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Ag^-crystal B 2 ' s if n = 4. Therefore, by the main results in [To| (resp. [lij when 
n — 4) the scattering rule is the same as for the (resp. A^) case, namely: 

Ti{ P ) = a A \ P ), 

for sufficiently large t. 

Assume the statement is true for < y 2 . By corollary 1, T^{p) = z kl ((si,0),b) ® 
z k 2 -i^y^ _|_ \^y 2 — so the inductive assumption holds. Therefore, T*(T^(p)) 

is a 2-soliton state, for sufficiently large t, and 

T t r {T^p))=H A \T^p)). 

Therefore, by Lemmas [3] and HI we have 

Also, by Lemma [3j we have b(lZ AS (p)) = b(p), and by Corollary [2] b (p) = [If] since 
p is a 2-soliton state and ?/2 > 0. By Lemma HI b(T r (p)) ® u r = TZ(u r eg) b(p)) = 
1Z(u r ® |~2~~|) — [2] Cg> M r . Therefore b{T r (p)) = [2] and, by repeated application 

of Lemma El we have 6(2? (p)) = [2]. Therefore, b{K AS {p)) = b(T l r (p)) = [2], 

and Tt|(T r *(p)) = 2^(72. Therefore, the crystal maps in (I2T1) can be inverted 

to yield 

T$(p) = TZ A \p) 

is a 2-soliton state. 

□ 



B. Examples of scattering of D, t -solitons 

Here we give some examples of D„ soliton scattering, for different n. 
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Example: n = 4 : 

t _ Q 111112211111111111111111111 
344224322222222222222222222 



t= 1 : 
t = 2 : 
t = 3 : 
t = 4 : 
i = 5 : 
t = 6 : 
t = 7: 



111111122111111111111111111 
222344243222222222222222222 

111111111221111111111111111 
222222344432222222222222222 

111111111114211111111111111 
222222222344322222222222222 

111111111111142111111111111 
222222222222333322222222222 

111111111111111122111111111 
222222222222223433322222222 

111111111111111111122111111 
222222222222222234233322222 

111111111111111111111122111 
222222222222222222342233322 



The initial state corresponds to z°((3, 0), (0, 3), (2, 1)) <g>2~ 5 (0, 2), (2, 0), (1, 1)). We 
compute the i?-matrix of (A^ } )® 3 -crystals as follows: ft aff (z°((3, 0), (0, 3), (2, 1)) <g> 
z- 5 (0, 2), (2, 0), (1, 1))) = z-\(2, 0), (1, 1), (0, 2)) ® ^((l, 2), (2, 1), (2, 1)). Ac- 
cording to Theorem [7] we expect this to correspond to the final state, which is the 
case. 

Example: n = 5 : 
t = 0: 

t = 1 : 



t = 2: 
t = 3: 
t = 4 : 
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1 


1 


1 


1 


2 


1 
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1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


i 


1 


1 


1 


3 5 


4 


3 


2 


2 


4 


5 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 
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2 


2 
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1 1 
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5 
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4 


5 


3 


2 


2 


2 


2 


2 


2 


the initial state corresponds to z° I 


(2,2), ( 


'2 

v0 


1 
1 


1 

2 


0* 


r 


)« 


5 2 


-6 


( 



(1,1), 
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j . We compute the .R-matrix for © A^-crystals as follows: 

^(( 2 , 2 ),e^3)-^(M),(j;-)))=^((M),(j;-)} 

z^ 1 ((2, 2), | g ^ ' wri ich again confirms Theorem [3 
Example: n = 6 : 

t = 
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t = 3 : 
t = 4 : 



In this example, the initial state corresponds to z ((3, 2), (0, 1, 1, 1, 0, 1, 0, 1)) © 
2T 7 ((0, 2), (0, 0, 0, 0, 0, 2, 0, 0)). We compute the combinatorial .R-matrix for A^ © 
/^-crystals as follows: ft Aff (z°((3, 2), (0, 1, 1, 1, 0, 1, 0, l))©z~ 7 ((0, 2), (0, 0, 0, 0, 0, 
2, 0, 0))) = z~\{2, 0), (0, 0, 0, 1, 0, 1, 0, 0)) © z°((l, 4), (0, 2, 0, 0, 0, 1, 1, 1)). This is 
in agreement with Theorem [JJ 
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